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Abstract 



We represent an integration algorithm combining the characteristics method and Hopf- 
Cole transformation. This algorithm allows one to partially integrate a large class of 
multidimensional systems of nonlinear Partial Differential Equations (PDEs). A specific 
■ generalization of the equation describing the dynamics of two-dimensional viscous fluid 

\Q \ and a generalization of the Korteweg-de Vries equation are examples of such systems. The 

1 richness of available solution space for derived nonlinear PDEs is discussed. 

1 Introduction 

■ 

A possible way to construct new classes of (partially) integrable nonlinear multidimensional 
Partial Differential Equations (PDEs) is combining features of different classical integration 
y—i ■ algorithms. Thus, the multidimensional nonlinear PDEs having features of PDEs linearizable 
(C-integrable pQ) by the Hopf-Cole transformation [21 E] and integrable by the Inverse Spectral 
Transform Method (ISTM) [H El El E] (S- integrable) have been derived in refs.JEl IS]- A com- 
bination of the characteristics method [10] with the Hopf-Cole transformation is proposed in 
refs.JTTJ IS]. The characteristics method was also combined with the ISTM in refs. [121 IS] and 
with the commuting vector fields [TBI IT4] in ref . [15] . There are also several generalizations of 
the characteristics method [161 El HH1 HS1 [20] integrating certain classes of systems of multi- 
dimensional nonlinear PDEs. In addition, there is a relation [21] between the nonlinear PDEs 
.^J [ linearizable by the Hopf-Cole transformation and integrable by the ISTM. 
^ ; The algorithm proposed in this paper allows one to partially integrate a new class of multi- 

dimensional nonlinear PDEs using the features of the characteristics method and Hopf-Cole 
transformation. Unlike the nonlinear PDEs derived in refs. [Ill IS], a new type of nonlinear 
PDEs involves additional arbitrary functions of fields and their derivatives. Note that the term 
"partially integrable PDEs" has two different meanings. First, the PDE is partially integrable 
if the available solution space is not full [22] • Second, the PDE is partially integrable if its 
solution space is described in terms of the lower dimensional PDEs |24j. Our nonlinear PDEs 
are partially integrable in both senses. 

The general result is covered by the following theorem. 
Theorem. Let W(x,t) be a solution to the algebraic system 

W(x, t)F(u(x, t),x, t) = 7](u(x, t)), (1) 

W(x,t)F u (n)(u(x,t),x,t) — r] u (n)(u(x,t)) = ri^ n \x,t), n—l,...,K, (2) 

det W ^ 0, (3) 

t = {t 1 ,...,t M }, x = {xi, . . . ,x N }, u = . . . ,u (K) }, (4) 
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where W is the N% x Ni matrix function of independent variables x and t; F, rj and rf n ' are 
the N\ x N 2 (Ni > N 2 ) matrix functions of arguments with 

Ni = N 2 (K + l). (5) 

Let rj and rf n > be arbitrary functions, while F satisfy the following system of compatible linear 
PDEs: 

M 

F tm - F x V^\u) = FTH( U ), m = 1, . . . , M, (6) 
j'=i 

where y( m -?) and r*"^, are arbitrary iV 2 x N 2 diagonal matrix functions of u. Then elements of 
the matrix function W satisfy the following system of nonlinear equations: 

N 

W t JVS[u) -^W^W-^V^iu) + (7) 
3=1 

K N 

E ( u t - £^^(«)y°0M) + ^) r(m) n = 0' 
t=i i=i 

where it*- 1 ) are arbitrary functions of x and t. 

Proof. First, we differentiate eq.flT]) with respect to t m and xf 

K K 



E<M := W tm F + WF tm + W^F^ufl = £ 7^)42, (») 

i=l i=l 

:= W Xn F + WF Xn + ^««S = E ( 9 ) 



8=1 



Consider the combination — Xljli E^^V^ 171 ^ and take into account eqs.(j2]) and ([6]). We 

obtain: 



N K N 



W tm F - J2 W X .FV^ + (42 - E u S^ (mi) ) + WFY^ m \u) = 0. (10) 

j=l i=l j=l 

Since is invertible, from eq.([T]) we have: F = W~ x r\. Substituting it into eq. fllOl) we obtain 
eq.©. □ 

Remark 1. If W = 1, rj = 0, JVi = 1, N 2 = K and det{(F u ) uU) : i, j = 1, . . . K} ± 0, then 
eq.([2D may be disregarded. Herewith eq.flTJ) yields the usual characteristics method [IU] and 
eq.(J7]) must be replaced with the system 



(») 

<m 



N 

-£<V<"*> = 0, i = l,...,K, m = l,2,.... (11) 

3=1 



If N\ = N 2 , rj = F Xl and det F 7^ 0, then W = F Xl F~ l so that we have the matrix Hopf-Cole 
transformation. To derive the nonlinear PDE for the matrix W in this case we have to replace 
eqs.Q by another system for F. For instance, if 

F tm = d^F + FB {m) - B {m) F, m = 1, 2, . . . , (12) 
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(where £?( m ) are diagonal constant matrices) then W satisfies the matrix Burgers hierarchy 

Remark 2. If N 2 = 1, i.e. K = N% — 1, then eqs.flT]) and (T5]) may be written as the following 
single matrix equation: 

WF = j), (13) 

where F and fj are block-row matrices 

F = [F F uW . . . F u (Ni-i)], fj = [rj 7/ 1 ) + VuW . . . ry^ 1 ) + 7^-4 (14) 
To satisfy the condition (|3]) we require det F^O and det f) 7^ 0. Then 

V^ = r)F" 1 . (15) 

Emphasize, that r/(u), r)( l >(x,t) and V^ m ^(u) are arbitrary functions of arguments with 
appropriate matrix dimensions. 

Note that functions vS 1 ' and rf*> (i = 1, . . . , X) are arbitrary functions of x and £ in eq.([71). 
Consequently, we have to add relations among vP\ rf % ' and elements of W. Further, each matrix 
equation with fixed m in the system (JTj) involves N\N 2 scalar equations for Nf scalar fields 

= 1, . . . ,Ni). Thus, the complete system consists of several equations ([7]) with different 
m = 1, . . . , T, so that N\N 2 T = Nf. Consequently, the nonlinear system is iV + T-dimensional. 

2 Nonlinear PDEs corresponding to Ni = 2, N2 = 1 

Consider the case Ni = 2, N 2 = 1, T = 2 and denote w 1 - 1 ^ = u. Then the (iV + 2)-dimensional 
system of nonlinear PDEs (J7|) reads: 

v 



^ := A™' 1 \^W al ) tm - J2(W al ) Xj V} mj \u)j + (16) 
A^ V $(W,u) (u tm -^u Xj V} mj \u)j +A ( V< m) = 0, a = l,2, 



m = 1,2. 



where 



= W227&i(u) - W 12 7&l(u), ^ = - ^21^l(«), (17) 

A v = (W n W 22 - W 12 W 21 ). 
Let 77 be independent on w for the sake of simplicity, 

u = W n , v = W 12 , p = W 21 , q = W 22 . (18) 
Then, introducing the operators acting on an arbitrary function h(x, t) by the formula 

N 

L^\h) = h tm -Y J h Xj Vl mj \u), (19) 



we write eg. ( 1161) as 

L^\u) + GnL { - m \v) + G 12 T ( f > = 0, m = l,2, (20) 

L {m \p) + G 2l L^\q) + G 22 L {m \v) + G 23 rS m) = 0, m = 1,2 (21) 

"11 = Ft! i (22) 

A w r/n 

G12 = tt\ ) (23) 

A (v) v m + A (w;l) 

= (24) 

G22 = „, . (25) 

23 A( w ^ A {w;l)( A (w;l) + 



where 



= - vt/2!, A (u,;2) = wfci - P?n, = - pu. (27) 



Owing to the arbitrary functions rj!-l\x, t) (i = 1,2), coefficients G\ 2 and G 2 s may be arbitrary 
functions of x and t. Consequently, we may impose the following relations 

G 12 = h(U), G 23 = f 2 (U) (28) 

where list of arguments U in arbitrary functions /j (i = 1,2) involves the fields u, u, p, g and 
their derivatives: [7 = {u, v,p, q, u x , v x ,p x , q x , u u v t ,pt, qt, ■ ■ ■ }• Substituting G\ 2 and G23 from 
eqs.(l23l) and f l26l) into eqs. ( 1281) and solving them for 77^, z = 1, 2, we obtain: 

■ff= ^-f' A . i = M. (29) 
Now eqs.([20]) and (|3TJ> read as 

L( m )( W ) + G 11 ^ {m) (^) + Ari m) = 0, m = l,2, (30) 
L^ m \p) + G 2l L^ m \q) + G 22 L^ m \v) + f 2 T { ; n) = 0, m=l,2, (31) 

where, in virtue of eqs. (l2Uj) . read 

Gu = Ag^ = A^-pr/n) 
AW7711 r]n(uq — pv) 
urj 21 - prj n 

(^21 



G 



22 



A^ 2 \A^f 2 -A^ V21 ) = (uri 21 - pyu)(vV2i(p - f 2 ) ~ gHgi - W2)) 

A^A^ 1 )/^! Vn(Q u ~ P v )(Q r ]U ~ v V2l) 



Eqs. (l30l) and (!30|) must be considered as the ultimate general (iV + 2)-dimensional system of 
nonlinear PDEs in this example. 
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2.1 Reduction of the system flSDl KID to the single PDE 

There is a particular reduction of the system (130] [31!) leading to the single nonlinear PDE for 
the function u. Namely, if the second term in the system (!30|) is negligible, then this system 
reads 

N 

«*» - E + A r ! m) = 0, m = 1, 2, (33) 

i.e. u satisfies the system of two commuting flows. In this case, it is reasonable to take fi as 
an arbitrary function of u and its derivatives, while eq.( l3ip may be disregarded. Of course, 
commutativity condition imposes additional differential relation on the field u. To avoid this 
relation we have to disregard eqs. (l3"3"|) with m — 2, keeping the only equation 

N 

«*-£^i (u) («) + /i I f ) = - (34) 

3=1 

This reduction may be realized by means of the multi-scale expansion with 

u ~v~p~q~eM, ^ mi) ~e (y) , d Xk ~ £ ( x \ <9 4m ~e (x) £ (y) , h ~ e w e (l) e w , (35) 

where all e's are small positive parameters (<C 1). In this case Gn ~ rM^. Consequently, 
the second term in eq. (130]) is of the order e^ u '(e^e^ v ') 2 , while all other terms are of the order 
£ (w) £ (x) £ (V) _ Thus, the leading e^'e^'e^ -order terms yield eq. (13"3"]) . 

Example 1: three-dimensional viscous fluid flow. Let N = 3, V^(u) = u, f = 

(A = Ylt=i &Xi i s the three-dimensional Laplasian). Then eq. (jM|) yields a special case of three- 
dimensional viscous fluid flow with the constant kinematic viscosity v 

3 

u tl — u Xj u — uAu = 0. (36) 

3=1 

A two-dimensional (M = 2) version of this equation was derived in [23] using a different 
algorithm. 

Example 2: multidimensional generalization of the Korteweg-de Vries equation 
(KdV). Let = u, V™ = (j > 1), / = ^=1^32^ ttijk = const. Then eq.® 

yields an (N + l)-dimensional version of KdV: 

TV 

U tl -U Xl U+ 22 a ijkUx z XjX k = 0. (37) 
i,j,k=l 

2.2 Solution space to the nonlinear system (I30U31I ) 

To obtain solutions to the derived system of nonlinear PDEs (I30|l3ip we, first, write expressions 
for the elements of the matrix function W given by the system (j!5p with N% = 2: 

W-i = V * F21 - ^ F ^ u w- 2 = v ii Fu ( 38 ) 
F , 2i(Fii) u — Fii(F 2 i) u ' F 2 i(Fu) u — F n (F 2 i) u 
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where F a \ (a = 1, 2) are solutions to eq.(E]): 

(M M \ 

*1 + tmV}^ (U),...,X N + J2 tmV^ («) ■ (39) 

m=l m=l / 

Substituting from eqs. (!38|) and 7$ (i = 1,2) from eq. (l29l) into the system (fl8|) we obtain 

VuFh. - fiF n F 2 i - fimi(F 21 ) u 
U fi(F 21 (F n )u-F 11 (F 21 ) u ) ' 1 > 

V21F2I — /2-Pll-^21 — flV2l(F2l) u . . 

P /i(F2i(Fii) u -F n (F 21 ) u ) ' 1 j 



-VuFn + fiF^ + fmuiFi 



l u 



fi(F 21 (F 11 ) u -F 11 (F 21 ) u ) ■ 
-V2iFu + f 2 Ff 1 + fiy 2 i(Fn)u 
fi(F 21 (F u ) u -F u (F 21 ) u ) ■ 



(42) 
(43) 



If fi depend on all fields and their derivatives, then the system (140]) is a lower dimensional 
system of PDEs for the fields u,v,p,q (dimensionality of this PDE must be less then the 
dimensionality of the system ( 130|3T|) ). Derivatives in eqs.( |40H43|) appear due to the functions fi 
(i = 1,2), which may depend on the partial derivatives of the fields u, v, p and q. In particular, 
if fi (i = 1, 2) are the function of single field u and its derivatives, then only eq. (14*0l is a PDE 
for the function u, while eqs. (14111431) are the non-differential equations for the fields v, p and q. 

Now we briefly characterize the solutions u, v, p and q to the system f l30f3~T|) . which are 
implicitly given by eqs. (1401143]) . As was noted above, the system (1401143 p is a system of PDEs for 
the fields u, v, p, q whose dimensionality is less then N + 2, i.e. we do not completely integrate 
the original system of PDEs (I30|3ip . but reduce its dimensionality. Further, eqs. (1401143]) involve 
two arbitrary functions J-n, i = 1,2, of jV variables. For this reason, the available solution 
space is not full (for the fullness, the solution space to the (N + 2)-dimensional system of four 
PDEs must involve four arbitrary functions of iV + 1 variables). Thus the system ( 130|3T]) is 
partially integrable by our algorithm. Since all fields are given implicitly by eqs.( )40ll43]) . they 
describe the wave breaking in general. 



3 Conclusions 

The algorithm proposed in this paper allows one to solve a new class of multi-dimensional 
nonlinear PDEs using the features of the characteristics method and Hopf-Cole transformation. 
The principal novelty of this algorithm is the presence of arbitrary functions of fields in eq.([7]), 
i.e. functions rj,rj^ and Thus, this algorithm possesses new features which do not 

appear in both method of characteristics and method of direct linearization by the Hopf-Cole 
transformation. The derived nonlinear PDEs are partially integrable since, first, one has to 
solve a lower dimensional PDE in order to construct solutions and, second, the available solution 
space is not full. The important question is whether the multi-scale expansion f ]33]) may be 
properly implemented in the proposed algorithm so that eq. (l34p would be exactly solvable. 

It is important that, among the derived nonlinear PDEs, there are such equations that may 
be considered as multidimensional generalizations of known integrable models. The physical 
applications of these generalizations become obvious provided that some terms are negligible, 
see eq.( 13~4"]) reduced from the system ( 13"U]I3T]) . 
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